Energy density and localization of particles 
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We consider the use of the energy density for describing a localization of relativistic particles. 
This method is consistent with the causality requirements. The related positive operator valued 
measure is presented. The probability distributions for one particle states are given explicitly. 
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The concept of particle localization has been raising 
intriguing questions almost from the onset of the rela- 
tivistic quantum theory. On the one hand, while there 
is a broad consensus on what the localization formalism 
should be and how it can be used, the formalism itself 
is still incomplete. Different authors propose different 
quantities to describe localization. On the other hand, 
detection events are the experimental basis of particle 
physics, a notion of local observables is a cornerstone of 
the algebraic field theory |Q of and the coincidence anal- 
ysis plays an important role in the general quantum field 
theory jj], ^| and quantum optics ||. Moreover, the very 
concept of particles takes somewhat a 'nebulous charac- 
ter' yj in the field theory in a curved spacetime. A way 
to substantiate the notion of particles and to reconcile it 
with the local framework of quantum field theory lies in 
the analysis of the (model) detectors' excitations [Q, |). 
Additional motivation for deeper understanding of local- 
ization comes from the relativistic quantum information 
theory [[|. Usually information is encoded into discrete 
degrees of freedom, but its processing and retrieval are 
performed in certain places in space and at certain in- 
tervals in time. They are localized, so there should be a 
way to describe this localization. 

The origins of the problem can be traced to the fact 
that not only is there no spacetime localization oper- 
ator (tJ, but also no unique spatial position operator. 
Moreover, the spatial position operator does not exist 
for massless particles with the spin larger than ^. Even 
when it does exist, there may be no probability current, 
or causality problems are present fj[ [s[ p[ [Tof. 

In this Letter we investigate the possibility of parti- 
cle's localization by means of its energy density, taking 
as an example a free scalar boson field. We give an ex- 
plicit construction of the positive operator valued mea- 
sure (POVM) that describes a localization and derive 
resulting probability distributions. The limitation of this 
approach are discussed in the last part of the Letter. Wc 
work in the Minkowski spacetime, and set g 00 = +1 and 

h = c = i. 

It is known that particles cannot be confined into a 
volume with a typical dimension smaller than 
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where (E) is particle's expected energy Q. Therefore, 
what is actually used is not an operator and its general- 
ized eigenstates, but a projection valued measure that is 
associated with it. 

At least some of the localization problems can be 
solved if 'observables' and their spectral measures are 
abandoned for the generalized measurements that are 
described by POVMs @ |§. A POVM constitutes a 
non-orthogonal decomposition of the identity by means of 
positive operators [jl3|. It is the most general method to 
associate probability distributions with density operators 
in quantum theory. Nevertheless, some open questions 
still remain and it is not easy to derive explicit expres- 
sions from the general principles of covariance, causality, 
etc. i|. 

Problems with position operators and POVMs have 
led to the construction of various alternative densities 
and currents. One of the popular approaches is to use 
the energy density, i.e., an expectation value of the 
renormalized stress-energy tensor T 00 (x,t) = H(x,t). 
Here T^ lL , is a stress-energy tensor, the Hamiltonian is 
H = J 7id 3 x and in the Minkowski spacetime without 
additional boundary conditions the renormalization is 
achieved by normal ordering. 

Consider a one-particle state 
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where p° = E(p) = yjm 2 +p 2 , (p|q) = (27r) 3 (2p°)<S(p - 
q), and J dn{p)\ij}\ 2 = 1. The inner product of two states 
is calculated as the inner product of their wave functions 
in the momentum representation, 



(3) 



A wave function in the configuration space is defined by 
a generalized Fourier transform of ip(p) as 



tp(x,t) = / d/j,(p)ip(p)e 



i(p-x— Et) 



(4) 



It is a well-known textbook fact that only momentum 
space wave function has a probability interpretation. A 
construction of Newton and Wigner Q , which introduces 
a new wave function |Q 
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is only a partial remedy. It could not be incorporated 
into a continuity equation, |l0|, If \ip NW (-x., t)\ 2 is 

taken to be a spatial probability distribution then there 
are states that violate causality pl| in the spirit of the 
Hegerfeldt's theorem that will be described below ||. 

Energy density is free from these deficiencies. For a 
state l^), it is 



E(x,t) = <*|r 00 (x,i)|vK) 
= |V^(x,t)| 2 + |^(x,*)| 2 



-m 2 |^(x,i)| 2 . (6) 



Lorentz transformation properties are built into this 
quantity by its definition, and it is positive. Let us con- 
sider possible causality violations. The Hegerfeldt's the- 
orem in its strongest version proves a superluminal speed 
for an exponentially localized particle. If the probability 
of finding it outside a sphere of radius R is bounded by 



Prob^ fl < C 2 exp(-2 7 i?), 



(7) 



where C is some constant and 7 > m, then the state 
will spread faster than light. However, it was shown by 
Barat and Kimball that no physical state can satisfy 
this bound. If E(x,t) satisfies it, then both \tp(x, t)\ and 
\d t ip(x, t)\ are bounded by Cexp(— 7i?). It implies that 
both V'(p) an d ip(p)/ E(p) are analytic functions in the 
strip of the complex plane that is bounded by at least 
|Im(p)| < to, which is inconsistent with the branch cuts 
in E(p) at |p| = ±im. Therefore, the energy density 
cannot be 'localized' enough to violate causality. 

Probability to find a particle in the spatial region A 
that is centered at the point x may be considered a nor- 
malized energy density 



P A (x,<) 



E(pc, t) 
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This quantity still cannot be taken as a probability, be- 
cause it is inconsistent with linearity of the quantum the- 
ory. Consider a state p which is a mixture of two pure 
states and |\J/ 2 ), 



p = api + (1 - a)p 2 



(9) 



According to Eq. (||) probability densities for the states 
Pi = are 

P*( x ' *) = /,t, 1 u 1 iTr x = ITTu ■ ( 10 ) 



tr pi H 



It is well-known that the probability domain is convex 
(see, e.g., [§, [[7)). Therefore, 

p p (x,<) = api(x,t) + (1 - a)p 2 (x,i). (11) 

On the other hand, 

tr p7i(x, t) 



P P (x,i) 



tr pH 

atr piH(x, t) + (1 — a)tr p 2 H(x, t) 
atx p\H + (1 — a)tr p 2 H 



,(12) 



and these two expressions are generally different. 

We can trace the difficulty to the fact that the most 
general way to obtain probabilities from density opera- 
tors is by the means of POVM, where a normalization of 
probability follows from the normalization of a measure. 
A 'manual' normalization, as in Eq. (||) brings the non- 
linearity that is revealed in Eq. (|l2|). However, it is easy 
to find a remedy. Since Hamiltonian densities commute 
for spacelike separated points, 



[H(x),H(x')} = 0, (a;-a;') 2 <0 



(13) 



the Hamiltonian commutes with the Hamiltonian density 
at each point. Therefore, we can construct a POVM 
element as 



A(x,i) = ff- 1 / 2 W(x,t)ff- 1 / 2 . 



(14) 



This formal expression glosses over a number of technical 
points. First, H^ 1 is ill-defined. However, its action 
is well defined when we restrict the Hamiltonian to the 
non-vacuum states. We also face a serious problem in 
that while 7Y(x, t) is positive when restricted to the one- 
particle states, it is not generally so. We address this 
issue at the end of this Letter. 

Calculation (W|A(x, t)\^) leads to the result which is 
similar to the energy density from Eq. ([|) 

p(x,t) = (*|A(x,i)|*> 

= |VV3(x,t)| 2 + |9 t ^(x,i)| 2 + TO 2 |^(x,i)| 2 , (15) 
where an additional E~ x / 2 factor is added to "0(p), 

VKp) 



and 



ip(x,t) = / dp(p)ijj(p)e 



i(p-x-Bt) 



(16) 



(17) 



Working in the momentum space representation, it 
easy to see that indeed 



J p(x,t)d 3 yi = 1. 



(18) 



The arguments of Barat and Kimball for the energy den- 
sity are equally well applied to our probability distribu- 
tion p, so Hegerfeldt's theorem does not present a para- 
dox. Operators A(x, t) have another nice trait: local 
commutativity, which is a desired property for a localiza- 
tion POV measures [Q. Since we work at the fixed time, 
two spatially separated regions Ai and A 2 , Ai n A 2 = 
are mutually spacelike. We define a POVM that gives 
a probability for a particle to be localized at the spatial 
region Ai as 



(19) 
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Then Eq. @ ensures that 

[A(A u t),A(A 2 ,t)]=0 



(20) 



It is interesting to note that while ip(p) = 
2~ 1 / 2 E(p)~ 1 ip NW (p), its position representation is used 
in a different manner. Moreover, for the states with a 
well-defined energy, AE/E <C 1, the 'correct' result of 
Eq. (^J) and the 'wrong' one of Eq. (||), augmented with 
the prescription of Eq. (|ll|) to handle mixed states give 
nearly identical probability distributions. 

The POVM of Eq. @ satisfies the standard list of 
requirements for localization operators 1 12 and solves a 
problem of localization in one-particle theory. However, 
the notion of localization that is based on the energy 
density cannot have a universal validity. There is no 
probability density 'in general', but only a probability 
density that is related to a specific detection scheme. 

A simple example when it runs afool is the Unruh effect 
j|, H [l8| . An accelerated detector that moves in the 
Minkowski vacuum responds as an inertial detector would 
if immersed into a thermal bath of the temperature 



T = a/{2Trk B ), 



(21) 



where ks is Boltzmann's constant and o the proper ac- 
celeration. A perturbative derivation of this result is sup- 
ported by the Green's function analysis. Green functions 
for the accelerated detector are identical with the ther- 
mal Green functions of the inertial one at the temper- 
ature T. However, the expectation of the renormalized 
stress-energy tensor is zero in both inertial and acceler- 
ated frames. 

In a more complicated settings a question of positivity 
becomes acute. Classical energy density is always pos- 
itive, which is to say that the stress-energy tensor for 
scalar field satisfies the weak energy condition (WEC) 
T^ u u^u 1 ' > 0, where it' 1 is a causal vector. In the frame- 
work of general quantum field theory |l], ^] it is impos- 
sible. There are states |T) that violate WEC, namely 
(T|T Al „u' 1 u I '|T) < holds [pL9|| , where T^ v now is a renor- 
malized stress energy operator. For example, squeezed 
states of electromagnetic [[| or scalar field have negative 
energy densities . It is known that even if WEC 

is violated the average WEC 



T^u u dT > 0, 



(22) 



still holds when the integral is taken over the world line 
of a geodesic observer (inertial observer in the Minkowski 
spacetime) p2[ . There are also more stringent quantum 
inequalities that limit the amount of the WEC violation. 
Instead of infinite time interval they deal with a sampling 
that is described by a function with a typical width to 
plj . Behaviour of fields subjected to boundary condi- 
tions is more complicated, but similar constraints exist 



also in these cases To our ends we need the anal- 

ogous inequalities to hold for a spatial averaging. This 
is, however, impossible. A class of quantum states was 
constructed for a massless, minimally coupled free scalar 
field (superposition of the vacuum and multi-mode two- 
particle states). These states can produce an arbitrarily 
large amount of negative energy in a given finite region 
of space at a fixed time Q. In this and similar cases 
the spatial averaging over part of a constant time surface 
does not produce a positive quantity. Consequently, the 
probabilistic interpretation of the energy density fails. 
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